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Biaxial strain in coherent epitaxial thin films can have a pronounced effect on the point defect 
profile in the film material. Detailed fundamental knowledge of the interaction of strain with point 
defects is crucial in understanding the stoichiometry and resulting properties of strained thin films. 
Here we investigate the effect of biaxial strain on the formation energy of cation and anion vacancies 
using MnO as a model system. Our density functional theory calculations show that, as expected 
from local volume arguments, compressive strain favours the formation of cation vacancies. Inter¬ 
estingly, we find that small compressive and tensile strains lead to ordering of the resulting holes 
along the in-plane and normal direction respectively, which should manifest in different anisotropic 
properties in the two strain states. 


I. INTRODUCTION 

Point defects in ionic materials are ubiquitous and of¬ 
ten strongly affect or even determine the material’s prop¬ 
erties. They can be classified as extrinsic, such as impu¬ 
rities or dopants, or as intrinsic, when the crystal lattice 
is locally imperfect for example by missing, displaced 
or interchanged atoms. Intrinsic defects occur even in 
the purest of materials as their formation is thermally 
activated. The relative concentration of intrinsic point 
defects is thus a function of the Arrhenius exponential 
factor of their formation energy AGform, the exact form 
depending on the respective defect equilibrium^. As the 
number of atoms and electrons may change during defect 
formation, the formation energy depends on the chemical 
potentials of the atoms and electrons in their respective 
reservoirs^ in addition to the internal energy cost to form 
the defect. This leads to the well-established effect of the 
environment on the defect profile of a material via the 
partial pressure of the constituent elements (see e.g. Ref 

3 ). 

We have recently shown using first principles density 
functional calculations that besides a reduced oxygen 
partial pressure in the surrounding environment, oxygen 
vacancy formation in manganite perovskites is also facili¬ 
tated by tensile strain, which can, for example, be applied 
using coherent heteroepitaxy in a thin film^. Enhanced 
oxygen deficiency under tensile strain was indeed exper¬ 
imentally observed in films of Lao.TCao.sMnOs^ as well 
as in cobaltite perovskites^. The effect of strain on the 
formation of oxygen vacancies can be seen as the con¬ 
verse of the concept of chemical expansion^: Depending 
on the charge of a point defect relative to the ideal lat¬ 
tice site, its formation is accompanied by oxidation or 
reduction of neighbouring transition metal ions, which 
lead to changes in ionic radius and manifest macroscop- 
ically as changes in lattice parameters. Conversely, if 
one imposes lattice parameters, for example by epitaxial 
strain, the locally expanded/contracted volume will fa¬ 
cilitate the reduction/oxidation of transition metal ions, 
leading to a larger concentration of point defects, corre¬ 
sponding to this changed oxidation state. We have fur¬ 
ther shown that epitaxial strain has the effect of order¬ 



FIG. 1. The cubic rocksalt unit cell of MnO. Arrows show 
the direction of magnetic moments on the Mn atoms. In the 
AFM-II phase moments of the same spin form (111) layers, 
two of which are shown with the green and blue planes. Color 
code: Mn=purple (large), 0=red (small). 


ing point defects on inequivalent sites in the lattice, with 
tensile strain strongly favouring oxygen vacancy sites for 
which the broken bond axis lies within the strained plane. 

In this work we study the response of the cation (man¬ 
ganese) vacancy profile to epitaxial strain, which is ex¬ 
pected to be favoured under compression. For complete¬ 
ness we also study the strain response of anion (oxy¬ 
gen) vacancies, which from volume arguments should be 
favoured under tensile strain. We choose the rock salt ox¬ 
ide MnO as our model system because its Mn^+ cations 
are readily oxidised and it is structurally simpler than 
related perovskite oxides. Our main findings are that, as 
expected, cation vacancies are favoured by compressive 
strain, with strain additionally inducing an ordering of 
the holes on sites around the defect. Against our expec¬ 
tation anion vacancies are not favoured by tensile strain 
but rather also by compressive strain, which stems from 
the formation of an F-center instead of charge localisa¬ 
tion on the transition metal site. 

Above the Neel temperature of 118 K MnO is a para¬ 
magnetic insulator with the ideal cubic rock-salt crystal 
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structure. Below this temperature the transition metal 
magnetic moments order in the antiferromagnetic (AFM) 
type II arrangement shown in Fig. 1 , which consists of 
ferromagnetic ( 111 ) planes, which align antiferromagnet- 
ically with each other^. The exchange interactions in 
this magnetic phase induce a magnetic stress^ leading 
to a compression along the [ 111 ] direction^^’^^, which 
manifests itself as a rhombohedral distortion of the unit 
cell. Experimentally reported distortion angles range 
from 0.43° to 0.62°^’i2d3 

The main point defects occurring in MnO are known 
to be Mn vacancies Using electrical conduc¬ 

tance measurements, the charge state of these Yun was 
determined to be the neutral vacancy {Y^^ in Kroger- 
Vink notation) at low oxygen partial pressure (P 02 ) 
low defect concentrations^"^. At intermediate P 02 the 
charge state changes to Y^^ and finally to Y^^^ at high 
P 02 - The predominance of Y-^^ at low P 02 was later 
confirmed using ionic transport measurements, but a de¬ 
fect cluster of four Y^^ and one Mn interstitial (Mn**) 
better fits the ionic transport data at high pos^^- ^ 
theoretical study of Vmii in MnO suggested that the O 
atoms around the vacancy acquire a magnetic moment 
and that the moments of the next nearest neighbour Mn 
ions are reduced compared to bulk-like Mn ions^^. These 
calculations were carried out in a 64-atom cubic supercell 
however, which, as will be discussed below, may lead to 
spurious results. Recently it was predicted that Yun at 
concentrations higher than 20% may render MnO a room 
temperature ferrimagnetic materiaU^. 


II. COMPUTATIONAL METHODS 

Our density functional theory (DFT) calculations 
were performed with the PBEsol exchange-correlation 
functionaU^ using the VASP code^^“^^. Wave functions 
were expanded in plane waves up to a kinetic energy 
cutoff of 550 eV. We used projector augmented wave 
(PAW) potentials^^’^"^ with Mn(3p, 3d, 4s) and 0(2s, 
2p) electrons in the valance. The rhombohedral distor¬ 
tion in MnO is severely overestimated by local (LDA) 
and semi-local (GGA) density functionals, the DET+U 
method^^ applied to the Mn d orbitals resulting in a 
much improved agreement with experiments^. Within 
this method, structural changes occur in the range U=0- 
2 eV whereas larger values of U do not further change 
the structure's. In our calculations we applied U=4.0 eV 
to the Mn 3d orbitals, which also results in correct redox 
energetics in MnO^^. Reciprocal space for the rhombo¬ 
hedral 4 atom (a,b,c=5.362 A, a, /d,7=34.258°) antifer¬ 
romagnetic (AEM-II) unit-cell was sampled using a F- 
centred 7x7x7 grid. Oxygen vacancy formation was stud¬ 
ied in a 64 atom (a,b,c=8.865 A, (a,/3,7=90.857°) super¬ 
cell where k-space was sampled using a 4x4x4 F-centred 
mesh. Cation vacancies were studied in a 512 atom 
(a,b,c=17.730 A, (a,/d,7=90.857°) supercell with only the 
F point used for reciprocal-space sampling. 


For strained film calculations we imposed the epitaxial 
constraint by setting one of the three cell angles to 90°, 
constraining the two in-plane lattice parameters and re¬ 
laxing the length and direction of the lattice vector nor¬ 
mal to the substrate at each value of imposed epitaxial 
strain. Internal coordinates in all calculations were re¬ 
laxed until forces converged below 0.001 eV /A. 

For cation vacancies multiple competing localised hole 
arrangement patterns exist. We achieved a systematic 
study of their relative energies by imposing different pat¬ 
terns via the occupation matrix of the Mn^+ atoms sur¬ 
rounding the vacancy^^. The lattice was first relaxed 
with these constraints to obtain the lattice distortion cor¬ 
responding to the desired hole order. The occupation 
constraints were then removed but the hole localisation 
pattern was retained due to the lattice distortion. 

Defect formation energies were computed according to 
equation 12 ^ 28 . 

AF/j = F/gtoi ^ ^ T T F/Qorr* (1) 

i 

Here Edef and Egtoi are the DFT total energies of the 
defective and stoichiometric supercells respectively at 
a given level of epitaxial strain. For charged defects 
(g % 0 ), a neutralising background charge was applied 
and the formation energy corrected by the Fermi energy 
(Efermi, the chemical potential of the electrons) and Ecorr, 
a correction derived using the sxdefectalign code^^’^^ to 
align the electrostatic potentials of the stoichiometric and 
charged defective cell. The Fermi energy is defined rela¬ 
tive to the valence-band edge and can take values from 
zero up to the band-gap of the material. 

In this work we denote defects and their charge states 
using Kroger-Vink notation^, which is of the following 
form: Mg . In this notation M refers to the defect species, 
which for the vacancies we consider here is V. The sub¬ 
script S indicates the lattice site the defect occupies and 
can be either a species of the perfect lattice or an intersti¬ 
tial denoted as i. The superscript C refers to the charge 
state of the defect relative to the lattice site it occupies. 
A charge of zero is indicated by x and positive and nega¬ 
tive elementary charges by • and ' respectively. Creating 
a Mn vacancy without adjusting the charge (q=0 in Eq. 
1 ) leads to a neutral vacancy on a Mn^+ site, which thus 
has a relative charge of -2e, which in Kroger-Vink nota¬ 
tion is Y-^^. This can also be regarded as the Mn vacancy 
introducing two holes into the material. Adding electrons 
to the system (q=l and q=2 in Eq. 1) will change the 
relative charge state of the vacancy to Y-^^ and re¬ 
spectively, which means that there is one hole or no hole 
respectively introduced into the material. 

During vacancy formation, rii atoms of a given species i 
are removed from the system, which is taken into account 
by their chemical potential The chemical potentials 
of Mn and O can vary within a range limited by the 
decomposition of MnO to metallic Mn {fiun = FMn, metal) 
on one side and the condensation of O 2 (/io = ^ 1 ^ 02 ) on 
the other side. The chemical potentials are related by 
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FIG. 2. a) Computed atomic structure showiug the relaxatiou 
of the O atoms arouud a Ymu vacaucy. The relaxatiou leads 
to a siguificaut shorteuiug of the Mu-0 bouds towards the 
Mu ious highlighted iu blue, b) Schematic molecular orbital 
diagram for the Mu-0 boud highlighted iu pauel a). The boud 
shorteuiug iucreases the euergy of the tt* molecular orbital, 
which becomes a prefereutial site for hole accommodatiou. 


MMn + Mo = MMnO,buik, and we express /iMn as a function 
of fiQ in this work. 


III. RESULTS DISCUSSION 
A. Mn vacancies in bulk MnO 

To analyse the electronic structure of a cation vacancy 
in MnO, we start with the simplest case: the 
cancy (which does not introduce holes in the lattice) in 
unstrained bulk MnO. In Fig. 2a) we show our calcu¬ 
lated two-dimensional section through such a The 

oxygen atoms surrounding the vacancy relax away from 
the vacancy, leading to shortened Mn-0 bonds between 
the relaxing O ions and six Mn ions, four of which are 
visible and highlighted in blue in Fig. 2a). This leads 
to an enhanced overlap and hence a higher energy of the 
TT* molecular orbitals formed between the t2g Mn 3d and 
the O 2p orbitals as schematically shown in Figure 2b). 
These states above the valence-band edge in the DOS 
have primarily Mn t 2 g character and are localised on the 
highlighted Mn ions. Due to their higher energy these 
states accommodate holes more readily than the corre¬ 
sponding molecular orbitals on regular Mn ions without 
shortened bonds. 

We find that larger Hubbard U values applied to the 
Mn d orbitals lower the energy of these occupied anti¬ 
bonding states, making their character gradually less Mn- 
like. We expect this to lead to a transition from well- 
localised polaronic holes to more delocalised band-like 
behaviour. We have determined however that unphysi- 
cally large values of U (> 8 eV) are required to induce 
a significant change in band-edge character, which im¬ 
plies that the picture of well localised holes is valid in 
Mn deficient MnO. 

This ordering of holes around the vacancy implies that 




FIG. 3. Mn vacancy formation energies (ambient condi¬ 
tions: iio = -4.85 eV, Efermi = 2.8 eV) in bulk MnO for dif¬ 
ferent charge states and arrangements of holes. Mn atoms 
accommodating holes are highlighted in blue. Note that the 
rhombohedral distortion as indicated by the angle a in panel 
a) is exaggerated for clarity. Ymu with two holes on a) two 
opposite Mn, b) two Mn in an angular arrangement with an 
angle > ^ and c) with angle < For Ymu fke single hole 
can be localised on d) one Mn or e) partially distributed over 
two opposite Mn while in f) for Ymu there is no hole. 


in order to fit even only one bulk-like Mn between Mn 
vacancies one needs a supercell with at least four cu¬ 
bic rock-salt MnO units along each of the crystal axes. 
This requires us to use a 512 atom 2x2x2 super cell of 
the cubic 64 atom AFM-II unit cell. Previous work us¬ 
ing smaller 64-atom cells^^ was likely affected by the fact 
that these hole-accomodating Mn ions along one direc¬ 
tion were equivalent by periodic boundary conditions. 
This prevents charge localisation as the hole count is arti¬ 
ficially doubled on a given Mn site. Indeed, if we perform 
computations in the smaller 64-atom cell, we recover the 
delocalised half-metallic state reported reported in Ref. 
16. 

Depending on the charge state of the cation vacancy, 
between zero and two holes are created. The state with 
two holes (Vj^^) leads to the lowest occupation of empty 
high-energy orbitals and is thus expected to be energet¬ 
ically most favourable compared to the y^n ^Mn 
with only one or no hole respectively. This is indeed what 
we see in Fig. 3 where we report the lowest formation 
energy for a y^n 0.415 eV (Fig. 3 a), whereas the 
formation energies for (Fig. 3 d) and (Fig. 3 
f) are much higher at 2.017 eV and 3.715 eV respectively. 
These formation energies are computed for a Fermi en- 
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ergy of 2.8 eV but the situation does not change through¬ 
out the whole accessible range of the Fermi energy. 

Since there are a maximum of two holes but six pos¬ 
sible Mn sites to accommodate them, different patterns 
of distributing the holes are possible. We also considered 
states where one hole was distributed over multiple of 
the six possible Mn sites, which we refer to as distributed 
arrangements. We attempted to stabilise all likely com¬ 
binations and were able to achieve the patterns shown 
in Fig. 3; other cases decayed into one of these pat¬ 
terns. In all cases the holes remained localised on the 
Mn sites highlighted in Fig. 2a) forming polaronic states 
with corresponding lattice distortion. While the under¬ 
coordinated O atoms relax away from the vacancy, short¬ 
ening the Mn-0 bonds from 2.21 A in the bulk to 2.11 
A as shown in Fig. 2a), the bonds which accommodate 
holes are additionally shortened to 1.92 A. 

For the y^n defect we see slight variations in the for¬ 
mation energy depending on the hole pattern. The most 
favourable hole arrangement pattern is when the two 
holes are on opposite sides of the vacancy (Fig. 3 a). 
Due to the slight rhombohedral distortion, two different 
angular arrangements of holes exist - one where the two 
hole-carrying Mn lie at an angle of more than 90 degrees 
with respect to the vacancy (Fig. 3 b) and one where 
this angle is smaller than 90 degrees (Fig. 3 c). We find 
that the situation with the larger angles is favoured by a 
minute amount, consistent with the electrostatic interac¬ 
tion between like-charge holes. Regarding the less stable 
defects, we note that hole distribution (half a hole 
on two opposite Mn) shown in Fig. 3 e) is less stable 
than the localised state (one hole on a single Mn) shown 
in Fig. 3 d). Hole distribution thus seems unfavourable, 
which is in agreement with the results for where all 
distributed patterns were found to be unstable and decay 
to fully localised orderings. 

From the above results on the formation energy of dif¬ 
ferently charged Vmii in bulk MnO it appears that the 
y^n das by far the lowest formation energy at all possi¬ 
ble values of the Fermi energy. This is in agreement with 
experimental findings at low P 02 ; tbe defect clusters re¬ 
ported to form at higher pos^^ were not considered in 
the present work. 


B. Mn vacancies in strained MnO films 

Having established the electronic structure of the Mn 
vacancy in bulk MnO, we will now consider strained thin 
films but restrict our analysis to defects, which, as 
established above, have significantly lower formation en¬ 
ergies than the other charge states. In a strained thin 
film, the likely hole arrangement patterns become more 
complex than in the bulk material, as strain lifts the sym¬ 
metry between the in-plane and out-of-plane Mn sites 
around the defect. Out of the resulting possibilities for 
hole ordering, we find that the partially distributed state 
with holes on two in-plane and two out-of-plane Mn sites 
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FIG. 4. Different likely arrangements of holes for Ymr (4=0) 
vacancy formation in a strained thin film: a) one hole on 
each of two opposite in-plane Mn, b) holes on two opposite 
out-of-plane Mn and c) partially distributed holes on four in¬ 
plane Mn. In panel d) the formation energies for these hole 
arrangements are shown. 

is unstable as is the distributed arrangement over six Mn 
sites. The remaining stable configurations are schemat¬ 
ically shown in Figures 4a) to c). From their respec¬ 
tive formation energies shown in Figure 4d) we see that 
the partially distributed hole arrangement on four Mn 
ions is less favourable than the two localised hole ar¬ 
rangements. The lowest energy hole arrangement is thus 
one hold on each of two Mn through the whole range of 
strains but switches from an in-plane pattern under com¬ 
pressive strain to an out-of-plane pattern under tensile 
strain. This transition could be measurable as a different 
anisotropy in the properties of Mn deficient MnO films 
under slight tensile and compressive strain. 


C. Oxygen vacancies 

Despite the fact that oxygen vacancies have limited ex¬ 
perimental relevance in MnO, we have studied the strain 
response of their formation energy for completeness and 
to see if it follows the local volume variation model^. 
Looking at the formation energy of the Vq* defects shown 
in Figure 5a) we see that 4% tensile strain very slightly 
favours oxygen vacancy formation. There is however a 
much larger lowering of the formation energy observed 
for compressive strain, which is in sharp contrast to the 
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FIG. 5. a) Oxygen vacancy formation energy (ambient con¬ 
ditions: iio = -4.85 eV) as a function of strain for the charge 
neutral (Vq*, q=0) oxygen vacancy. Electronic densities of 
states for b) stoichiometric MnO, c) a Vq* oxygen vacancy 
with two electrons forming a defect state 1.2 eV above the va¬ 
lence band edge. The DOS have been aligned at the valence 
band edge and the Fermi energy is indicated by the dashed 
line. The 0.004 e/A^ isosurface in panel d) shows the charge 
density associated with the F-center defect state for the Vq*- 


behaviour expected for defect electrons localised on tran¬ 
sition metal sites^. 

From the electronic density of states we see by com¬ 
paring the stoichiometric case (Fig. 5b) with that for 
Vq* (Fig. 5c), that oxygen vacancy formation leads to 
the appearance of a two-electron defect state about 1.2 


eV above the valence band edge. The isosurface of the 
charge density associated with this defect state (Fig. 5 
d) shows that the additional charge is located on the va¬ 
cancy site, forming an F-center. The formation of an 
F-center is not surprising as the reduction of Mn^+ to 
Mn^+ would be highly unfavourable. The markedly dif¬ 
ferent strain dependence of the oxygen vacancy formation 
energy can thus be attributed to the formation of this F- 
center state. A more detailed study of different charge 
states and charge localisation will be presented elsewhere. 


IV. CONCLUSIONS 

Using density functional calculations we have studied 
both Mn and O vacancies in epitaxially strained MnO 
thin films. We show that neutral y^n defects, which in¬ 
troduce two holes into the material, are significantly more 
favourable than charged defects in which these holes are 
filled with electrons. This is consistent with the fact that 
the additional electrons associated with charged defects 
occupy high-lying anti-bonding states. Strain has the 
expected effect on the formation energy of such a lattice 
contracting defect, making it more favourable under com¬ 
pression. Interestingly strain favours different arrange¬ 
ments of holes around the defect, which should lead to 
different anisotropic properties under slight compressive 
and tensile strains. Due to the formation of F-centers in¬ 
stead of charge localisation on the transition metals, oxy¬ 
gen vacancies are also favoured under compressive strain, 
in contrast to the simple chemical expansion model. 
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